Repeated injections of energy in the first 600 ms 
of the giant flare of SGR 1806-20 
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£> . The massive flare of 27 December 2004 from the soft 7-ray repeater SGR 1806- 

20, a possible magnetar 1-3 , saturated almost all 7-ray detectors 4-7 , meaning that 
the profile of the pulse was poorly characterized. An accurate profile is essential 
to determine physically what was happening at the source. Here we report the 
unsaturated 7-ray profile for the first 600 ms of the flare, with a time resolution 
of 5.48 ms. The peak of the profile (of the order of 10 7 photons cm -2 s -1 ) was 
reached ~ 50 ms after the onset of the flare, and was then followed by a gradual 
decrease with superposed oscillatory modulations possibly representing repeated 
energy injections with ~ 60 ms intervals. The implied total energy is comparable 
to the stored magnetic energy in a magnetar (~ 10 47 erg) based on the dipole 
magnetic field intensity (~ 10 15 G), suggesting either that the energy release 
mechanism was extremely efficient or that the interior magnetic field is much 
stronger than the external dipole field 2 . 
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At the onset of the giant flare of SGR 1806-20, plasma particle detectors on the Geotail 
spacecraft detected an extremely strong signal of soft 7-ray photon fluxes (integrated above 
~ 50 keV) during the initial intense phase of the giant flare (£=0-600 ms). Figure 1 shows the 
count profiles of two detectors, Amcp ( re d symbols) from the microchannel plates (MCPs), 
and Nqem (blue symbols) from channel electron multipliers (CEMs), where Acem is scaled 
by a factor of 280 to account for the sensitivity difference. The onset time (t—0) corresponded 
to 21 h 30 min 26.35 s Universal Time (UT), which was consistent with the expected arrival 
time of the onset signal at the Geotail position. To understand how the flare energy release 
occurred, the detailed time profile of the flare is important. Before the onset Amcp was at the 
background level (<~11 counts, shown by a black arrow), and then increased to 839 counts 
within 5.48 ms, so that the e-fold time of the initial rise was shorter than ~ 1.3 ms. After 
the intermediate level of 1,330 counts at £=5.48 ms, the MCPs were saturated and Amcp 
could not be determined until ~ 176 ms. Between £=22.7 and 170 ms, we could obtain Acem 
values instead. The scaled Acem increased to 25,900 at £=33.6 ms, thus giving an e-fold time 
of 9.5 ms. (A data gap between £=11.0 and 22.7 ms shown by the leftmost grey bar was due 
to the scheduled instrumental operation and not caused by the flare itself.) Between £=33.6 
and 55.5 ms the scaled N C em stayed at the peak level of ~ 25, 000-27, 000. After £=61 until 
170 ms the scaled Acem decreased gradually with oscillatory modulation, which suggests 
repeated energy injections at ~ 60 ms intervals. (Note that a similar injection profile was 
also seen during the impulsive phase of the giant flare of SGR0526-66 on 5 March 1979 
(ref.8).) After 176 ms Amcp became available again and showed a continuing exponential 
decay with an e-fold time of ~ 66 ms until £=380 ms. The decay profile of scaled A*cem 
available for £=210-308 ms is consistent with that of iV M cp. Between £=397 and 500 ms 
several humps were seen on the profile of Amcp- Although less significant, scaled Acem 
showed a similar hump for £=402-451 ms. (Note that the same humps were also detected by 
the BAT detector on the Swift spacecraft 7 .) The physical origin of these humps is not clear 
at the moment, but may represent some additional energy- releasing process. After £=470 
ms, A^mcp again decayed with an e-fold time of ~ 57 ms. 

To convert the observed count rates to physical quantities such as energy flux, we need 
the energy spectrum information, which was not available from the Geotail observation alone. 
We have therefore taken three reported function forms at the peak of the giant flare 4,5 ' 7 (Table 
1) and integrated them above 50 keV. We then found that the resultant estimations of photon 
number flux, energy flux and fluence (for £=0-600 ms) are almost independent of the choice 
of the energy spectrum, and are ~ 2.5 x 10 7 photons cm~ 2 s^ 1 , ~ 20 erg s^ 1 cm~ 2 and ~ 2 erg 
cm -2 , respectively. The corresponding total energy radiated from SGR 1806-20 is estimated 
to be ~ 5 x 10 46 (tt/Air) d 2 5 erg, where d i5 is the distance scaled by 15 kpc and f2 is the solid 
angle of the radiation. Here we note that the solid-angle factor (Q/Air) is not likely to be 
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as small as 1CT 2 , as is typically assumed for relativistic jets for GRBs. (il/4n)> 0.1 is more 
likely for the intense initial spikes of SGRs because they have been seen in all the three giant 
flares of SGRs that could have been detected without them. Therefore the presence of a very 
efficient mechanism is implied, which promptly releases (on a timescale of ~ 60 ms) a major 
fraction of the stored magnetic energy in a magnetar, E ma , g w 1.7 x 10 47 Bf 5 Rg erg (where 
B15 and R6 are the internal magnetic field scaled by 10 15 G and the radius of the magnetar 
scaled by 10 6 cm). Alternatively, as suggested by ref. 2, the internal magnetic field could 
be as strong as (5-10) xlO 15 G so as to permit the emission of multiple giant flares over the 
lifetime of a magnetar. 

As we noted above, there were humps in the light curve of the 2004 giant flare at 
£=400-500 ms. Similar humps were also observed 9 at 200-600 ms after the onset of another 
giant flare of SGR 1900+14 on 27 August 1998, whose total energy is smaller by a factor of 
~ 100 than the 2004 giant flare. From similarities in the timings of the humps despite the 
large difference in the total energies, we suggest that the humps more probably represented 
continuing energy injections, rather than the results of interactions of the flare ejecta with 
environmental matter. 



Methods 

The Low Energy Particle (LEP) experiment 10 aboard Geotail consists of an ion detector 
with seven independent MCPs and electronics systems, and an electron detector with seven 
independent CEMs and electronics systems, both of which are designed to measure plasma 
particles in the solar wind and magnetospheric environment. When the giant flare occurred, 
Geotail was at (-1.5997xl0 5 ,-97945, -19671) km using the Geocentric Solar Inertia (GCI) 
coordinates (J2000), which was in the solar wind about ~ 10 earth radii upstream from the 
bow shock. Although MCPs and CEMs kept measuring the solar wind ions and electrons 
throughout the giant flare interval, these particles were being selected electrostatically and 
came into the detectors mainly at some limited timings that fortunately did not overlap 
the giant flare interval: the contribution of solar wind ions to Amcp, which is the sum of 
counts over seven MCPs, was at most 30, and thus was negligible for the study of intense 
7-ray photons. On the other hand, during subintervals (t < 10 ms, 175 ms< t <200 ms, 320 
ms< t <400 ms and 450 ms< t <600 ms) , A*cem, which is the sum of counts over seven 
CEMs, was affected by solar wind electrons, and was not available for the 7-ray photon 
detection. 

Another fortunate factor was that the angular distance between the Sun and SGR 1806- 
20 was ~ 5 degrees, so that previous knowledge of the 'calibration' of the LEP detector to 
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be used as a soft 7-ray photon counter on the basis of solar flare photon analysis was directly 
applicable to the interpretation of the observed characteristics of photons from SGR 1806-20. 
By comparing the count rates of MCP and CEM with the hard X-ray 11 and 7-ray 12 data 
from the Yohkoh space solar observatory during major solar flares in 1997-2001 (ref.13), we 
have seen that MCP and CEM are sensitive to soft 7-ray photons above ~ 50 keV, where 
their sensitivities are evaluated as the product [e S] of quantum efficiency e and effective 
detection area S summed over seven MCPs and seven CEMs against 7-ray photons. (Here 
e is defined to include not only the detector response itself but also the attenuation factor 
inside the spacecraft.) From the spectral information provided by the Yohkoh observations, 
we have calculated photon fluxes in two energy ranges, L (50-500 keV) and H (above 500 
keV), and then estimated [eS]yicP,L and [€S]mcp,h separately for these two energy ranges. 
(This separation is possible because the energy spectra of incident solar 7-ray photons differ 
from event to event.) Figure 2 shows the calibration result. Along the line of sight towards 
SGR 1806-20, we have found that fe 5y MC p,L=(0.19±0.06) cm 2 and /e S/mcp,h=(0.22±0.16) 
cm 2 where systematic errors are included. The estimation of [eSjcEM has been also done 
with solar flare 7-ray photons, and the result is summarized as [e S , /cem~1/280 of /cS'/mcp- 
The smallness of [sSJcem as compared with [€S]mcp is consistent with the difference in 
physical sizes of CEMs (millimetres) and MCPs (several centimetres). 

It is noted that [eSjucp obtained above is by a factor 10~ 2 -10~ 3 smaller than those 
of conventional 7-ray detectors. Nonetheless MCPs were saturated during ~ 150 ms of the 
onset of the giant flare of SGR 1806-20. We have made that the dead-time analysis for 
MCPs and found that the characteristic dead time is ~ 4.3 yus, which is consistent with 
the pre-flight calibration of the LEP system as well as the calculated circuit time constant. 
Nmcp, shown in Fig. 1, is after the dead-time correction, which becomes significant above 
~ 1, 000 counts. On the other hand, CEMs, which are two orders of magnitude less sensitive 
than MCPs, were found to be free from the saturation effect even at the peak of the giant 
flare. 
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Table 1: Estimations of energy flux, fluence, luminosity and total energy for several models. 
(The Planck distribution with kT=175 keV is from ref. 4; the power-law distributions with 
exponential cut-offs at 480 keV and 800 keV are from refs. 7 and 5, respectively.) 
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Fig. 1. — Observed photon counts during the first 600 ms of the giant flare. iV MC p (red 
squares) and iVcEM (blue solid squares) show the counts of MCP and CEM instruments 
accumulated over a bin of 5.48 ms duration. Averages over three successive bins (16.44 
ms) are taken for A^cem (blue open squares) after 200 ms. Vertical error bars represent 1 a 
statistical deviations, which become smaller than the symbol size for Nmcp > 100. Grey bars 
show the data gaps every 32 bin (or every 187.09 ms) caused by the scheduled instrumental 
operation. Before t=0, the photon counts were below ~11. The black arrow shows this 
upper limit. The inset shows t=0-60 ms with a linear scale. 
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Fig. 2. — Calibration of the 7-ray sensitivity of the MCP detector. The observed MCP 
counts (counts per s) are compared with the counts synthesized from the solar 7-ray photon 
observations with the estimated [sSJmcp,!, and [€S]mcp,k for five solar 7-ray flares (Goes 
class X3.7 on 22 November 1998, X4.9 on 18 August 1998, X2.3 on 24 November 2000, X5.3 
on 25 August 2001 and X9.4 on 6 November 1997). The blue and red parts of each bar 
represent contributions from the photons in the energy ranges of 50-500 keV and above 500 
keV, respectively. 



